Inertial Josephson Relation for FIR Frequencies 
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Consideration of the balance of forces on superconducting condensate at low frequencies leads 
to the well-known Josephson Relation. Using the Ginzburg-Landau expression for the current, 
an expression relating the electric field to the vortex velocity via the magnetic field is obtained. 
This result is the Josephson Relation, supplemented by a term accounting for the inertia of charge 
carriers. This Inertial Josephson Relation may be used at all frequencies and may be viewed as 
the Josephson Relation extended to the case of sub-gap high-frequency response. When applied to 
vortex dynamics it yields the same conductivity as solution of the Ginzburg-Landau theory. 
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Vortices in the superconducting condensate are stable 
objects, which interact to form vortex matter [l|, 0]. The 
response of vortex matter may be studied at various fre- 
quencies. One potential practical application of vortex 
matter is as a memory device; one might imagine that 
data could be encoded in the positioning of vortices. To 
attain this degree of control over vortex matter, it must 
be clear theoretically how forces and fields interact with 
vortices at the desirable high frequencies of such a puta- 
tive memory device. 

Low-frequency vortex motion can be observed by mag- 
netic scanning [3[, but high frequencies are invisible via 
this method. Fortunately, in far infra-red (FIR) spec- 
troscopy it has been possible to make high-frequency ob- 
servations of the electric current to which the vortex ve- 
locity responds Although these data have then been 
interpreted using the Josephson Relation Q , one may ask 
whether this is legitimate at high frequencies. 

An alternative way to interpret the FIR data is based 
on the Time-Dependent Ginzburg-Landau (TDGL) the- 
ory 0] • These two approaches lead to contradictory 
results; there is a difference in phases which is illustrated 
schematically in Fig. Q] Actual experimental data is 
shown in Fig. [2j Both approaches can reproduce the ex- 
perimentally established frequency-dependent phase dif- 
ference between the electric field and the current, but 
they yield different vortex velocities. According to the 
Josephson Relation 

E = --v L xB (1) 

c 

the vortex velocity vl is in phase with the electric field 
E. In contrary, in the solution 0] of the TDGL equation 
derived near the transition line, the vortex velocity is in 
phase with the current. 
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FIG. 1: Phases of electric field E, current j and vortex velocity 
vl as functions of time t 

These different phases of the vortex velocity lead to 
different interpretations of the experimental data. The 
motion of vortices is given by the balance of forces (per 
unit length) 

-jxB-w 1 =i ) v L) (2) 
cB 

where the left-hand side includes the Lorentz force due 
to the mean current j and the pinning force proportional 
to the vortex displacement and the right hand side is 
a friction force. Identification of these forces corresponds 
to the Gittleman and Rosenblum (GR) model 0] which 
is sufficient for our discussion. Using the Josephson Rela- 
tion ([I]) one must interpret the experimentally observed 
phase difference between the field and current as a man- 
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ifestation of the pinning force, v ^ 0. Interpreting this 
same case using the TDGL theory, one would not require 
the addition of a pinning potential. 

The goal of the present paper is to show that the de- 
screpancy discussed above is actually a consequence of 
the inapplicability of the Josephson Relation at high fre- 
quencies. Using the current from the GL theory, we will 
show how one can derive the Inertial Josephson Relation, 
which is similar to the Josephson Relation except for a 
term important only at high frequencies, which is propor- 
tional to the acceleration of charge carriers. The Inertial 
Josephson Relation has in fact been known in the litera- 
ture for some time | L 8, 9] , but has been derived by much 
less direct means than that which we propose here (8| . In 
d,Q it is motivated by the application of hydrodynamics 
to the superfluid state fioj . 

Although our arguments are in many respects more 
widely applicable, for the sake of simplicity we re- 
strict our attention to a two-dimensional sample with 
perpendicular magnetic field B creating the triangular 
Abrikosov vortex lattice. We will derive the relation for 
an ideal sample with no pinning centers. Electrons are 
driven by a FIR light which we represent by a homoge- 
neous electric field parallel to the sample E = Ecoswi. 
The sample is of infinitesimal thickness, therefore cur- 
rents in the sample have negligible feedback effect on the 
electromagnetic fields acting on the sample. We neglect 
contributions important for the Faraday rotation or the 
Hall voltage. Finally, we work in the vector gauge with 
zero scalar potential, allowing us to avoid a distracting 
discussion of the electrostatic and electrochemical poten- 
tials. 

In principle one must solve for the order parameter ip 
from the TDGL equation 
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and then evaluate the current 



J = — 
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In the limiting case of linear response these two equations 
can be mapped to equations @ and (P) , respectively. As 
shown in 11[ , the TDGL equation © implies the balance 
of forces on vortices @ with v = 0. A second indepen- 
dent equation must be derived from the current formula 
(|3J). The electric field appears in the time derivative of 
the current, since E = —(l/c)(dA/dt). 

In the linear approximation the vortex lattice is undis- 
torted, and the order parameter reflects this as a func- 
tion of time. We denote by ipo (r) the order parameter in 
the equilibrium and by ip(v,t) its time-dependent value 
driven by the light. The coherent movement of the vor- 
tex lattice implies ip(r,t) = e _ ^ CL ' r - ) V'o( r — r L), where 
tl is a displacement of the vortex lattice at time t from 



its equilibrium position. The phase factor is given by the 
relative position of the displaced vortex lattice and the 
center of the vector gauge, C L = (7r/$ )B x r L . Here 
$o = 2nhc/e* is the elementary flux. When calculating 
03 / dt, the time derivative can thus be written in terms of 
space derivatives; d^p/dt — i(27r/$ )( v L • A)ip— (vl • V)^, 
where we have used vl = (dr^ / dt) and assumed vl par- 
allel to rL. We find 
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The transport current j is the mean value obtained 
by averaging over the elementary cell of the lattice. We 
thus average the time derivative of the current ([5]) and 
express the mean current via the Cooper pair velocity, 
j = e'nv, Here n = (|Vf) = (B/$ ) J ccll cfr|^| 2 is the 
mean density of Cooper pairs. 

e* at 2e*n x ' 
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Simplifying, 
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The total derivative of the current is zero under integra- 
tion, since the current itself is periodic on the lattice. 
Using the identity ba • c — ca • b = a x (b x c) we thus 
obtain the relation 



m * 9v s „ 1 „ 
^= E+ c VLXB 
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which is the Inertial Josephson Relation; we view it as 
extending the validity of the Josephson Relation (JXJ) into 
the FIR region. 

In Kopnin [ll| is given a derivation of the Josephson 
Relation ([T]) under certain rather broad assumptions. It 
is useful at this stage to make contact with this well- 
known result, and show how this is commensurate with 
the Inertial Josephson Relation <JSj> - Differences between 
11] and the present paper include Kopnin's use of gauge- 
invariant quantities, lack of assumption of homogeneous 
electric field, and use of the phase \ 01 the order pa- 
rameter ip = \ip\e lx explicitly in his equations. None of 
these differences is significant in the present context. We 
choose to make a certain choice of gauge, but the final 
result is one between physical quantities. The simplicity 
of our calculation is increased slightly due to our assump- 
tion of a homogeneous electric field, since we may take 
E outside (■ • • ) when calculating the averaged current 
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j. Kopnin uses the x field, which is subject to nonan- 
alytic behaviour due to a coordinate singularity when 
\tjj\ = 0, but this is done carefully and no descrepancy is 
to be found; we use the complex order parameter ip itself, 
which is not singular. 

The essential difference is as follows. Kopnin assumes 
the vortex configuration to move coherently at constant 
velocity, so that in the notation of [ll| (except that there 
A is used in place of tp), the vector potential is taken 
to be A(r,i) = A sta tic(r — vtJ) + Ai and the order pa- 
rameter ip(r,t) = V'staticlr — vl£) + Vh where Ai and 
■01 are taken to be small corrections. No compensating 
gauge transformation is needed here as the calculation 
is not performed in a particular gauge. In fact, in the 
time-dependent case these would need to be modified to 
read 



and 



A(r,t) = A 8ta tic(r-r L (t))+Ai (9) 



f/}(r,t) =-0 s tatic(r-r L (i))+0i (10) 



In addition, the longitudinal part of Ai may not be ne- 
glected, and we must insert the current 



j = -i(A 1 |^| 2 )^l 

c m* 
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With these modifications, relaxing the assumption of 
constant- velocity motion of the vortices, the derivation of 
Kopnin reproduces precisely the Inertial Josephson Re- 
lation (0. 

Let us illustrate in experimental context how the Iner- 
tial Josephson Relation ([5]) leads to conclusions different 
from those based on the Josephson relation (fTJ). In Fig. [2] 
is shown the imaginary part of conductivity observed 
with Fast FIR Transmission Spectroscopy by Ikebe et 
al With j = Re(je~ 4W *) the complex conductivity is 
defined by j = ctE. In the GR model, balancing forces 
as in ([2]) and using the Josephson Relation ((T|) to inter- 
pret the experimental data, one concludes that is =/= 
and there obtains an imaginary part to the conductiv- 
ity. From the TDGL point of view, again using ([2]) but 
with the Inertial Josephson Relation ©, one obtains a 
complex conductivity without recourse to postulating a 
pinning potential. 

While excellent experimental agreement has been at- 
tained through more empirical methods such as the in- 
terpolation between normal and superconducting states 
of Coffey and Clem [l| , we would like to stress that our 
derivation of the IJR relies only on the same basic phys- 
ical principles as the Josephson Relation. 

Using ?l = ivl/w, wc substitute the vortex velocity vl 
from ([2]) into the Josephson Relation ([T]), which yields 
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FIG. 2: Imaginary part of the conductivity for various mag- 
netic field strengths: The points are experimental data of 
Ikebe et al The Gittleman and Rosenblum model, us- 

ing the Josephson Relation, is plotted with pinning v — 5.9 
N/cm 2 . The result from TDGL theory is coincident with the 
use of the Inertial Josephson Relation, without pinning. 



Since the magnetic field is perpendicular to the current, 
the double vector product becomes E = $o-Bj/((?7 + 
iv/lu)c 2 ). The conductivity 



CTJR = 



(rj + iv/<jj)c 2 
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has non-zero imaginary part exclusively due to the pin- 
ning; v ^ 0. As one can see in Fig. [3J this model allows 
for a qualitatively good fit of experimental data since 
the imaginary part decreases as \/uj and increases with 
decreasing magnetic field. 

Now let us turn to the application of the Inertial 
Josephson Relation §E§ and show that the inertial term 
leads to complex conductivity even in absence of pin- 
ning. We evaluate the vortex velocity vl from the bal- 
ance equation ([2]) with v = and substitute it into the 
Inertial Josephson Relation ([8]). The resulting electric 
field is E = (<fr B/{r]c 2 ))j - «w(m*/(e* 2 n))j giving a con- 
ductivity of the Drude type 
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where B C 2 is the upper critical field at zero temperature 
and t = rim* c 2 / (e* 2 n$oB C 2) is the relaxation time of the 
condensate. 

In Fig. [2] we compare the imaginary part of conductiv- 
ity (|14j) with experimental data of Ikebe et al [IH using 
the relaxation time r = 7r 4 K 2 $oCN/(14C(3)c 2 i? C 2) in the 
dirty limit (llj . The experimentally established normal 
conductivity is o>j = 2 x 10 4 f2 _1 cm _1 and the upper 
critical field B c2 = 12. 2T, as in (l^. 

The GL parameter is n — 38, usual for 
NbN. From r we find the friction coefficient r/ = 
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7r 4 $ 5 C 2n<7N/(14C(3)noc 2 ). We have used the Lon- 
don penetration depth at zero temperature A 2 = 
m*c 2 /(27re* 2 n ), the GL coherence length £ = A/ft, and 
B C 2 = < i > o/(27r£ 2 )- The density of Cooper pairs we eval- 
uate from the TDGL solution in the small signal limit; 
n/n = (1 - T/T c - B / B c2 )/ (2/3 A ) with the Abrikosov 
coefficient (3a = 1.16 for the triangular lattice. Us- 
ing C(3) = 1.202 we obtain a = 2.495 cr N (l - T/T c - 
B / B C 2) / (B / B C 2 — iwr). A good agreement of the con- 
ductivity (fT4|) with experimental data shows that the set 
of equations @ and © captures the FIR response with- 
out any need for the addition of pinning to the theory. 

Let us revisit Josephson's derivation of relation (fTJ) 
in light of the present discussion. Josephson assumed 
a steady motion of vortices and condensate. Since the 
condensate is not permanently accelerated by passing 
vortices, he concluded that the electromotive force is 
balanced by a gradient of the electrochemical potential. 
Josephson's relation (fT]) represents the balance of these 
two forces. Away from low frequencies, it becomes im- 
portant that the condensate is continually accelerated by 
the alternating field. The balance of forces then leads to 
the Inertial Josephson Relation (|8]). Finally it is impor- 
tant to note, as did Josephson, that the E' appearing in 
his relation is not exactly the electric field E appearing 
in the present paper; the subtle difference is that Joseph- 
son's E' depends on the chemical potential rather than 
the electromagnetic scalar potential. 

The derivation here is restricted to the case of an in- 
finitesimally thin layer. It is reasonable to expect that 
its validity is more general, although the above analysis 
benefits from the simplicity of averaging over the uniform 
electric and magnetic field configuration. We did not con- 
sider the conductivity contribution of normal electrons; 
since the imaginary part of the normal conductivity is 
negligible, this is not important for our discussion. 

We have made use of Ginzburg-Landau theory to pro- 
vide a careful derivation of the Inertial Josephson Rela- 
tion, which relates electric and magnetic fields to not only 
vortex velocity, as in the Josephson Relation, but also to 
the acceleration of the charge carriers. We have shown 
that the same reasoning which leads to the Josephson 
Relation at low frequencies can be used to carefully de- 
rive the more generally applicable inertial version, first 



suggested by Abrikosov et al 0. 

We have illustrated the use of the Inertial Josephson 
Relation in interpreting experimental data in the con- 
text of high-frequency response. While the introduction 
of pinning is required to explain the data if the usual 
Josephson Relation is used, we have argued that in fact 
this use is erroneous; the Inertial Josephson Relation 
should be utilised at high frequencies, and this then does 
not require pinning. The analysis suggests that it is the 
inertial term which is responsible for the imaginary part 
of the conductivity in the high-frequency region. 
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